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Kepler's laws 




• Satellites orbit the earth following the same laws that govern the 
motion of the planets around the sun 

• Empirical laws descripting planetary motion 

• Kepler derived empirically three laws describing planetary motion 

• Kepler's laws apply quite generally to any two bodies in space which 
interact through gravitation (earths primary, satellite -> secondary) 
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Kepler's first law 



* "the path followed by a satellite around the Earth will be an ellipse" 



r-a-■i--a--I 

Flgu ne 2.1 The fori F 2 and F 2 r the semi major axis a, and the 
semiminor aids b of an ellipse. 

• The semi-major and semi-minor axes of the ellipse (a and b) are 
related to eccentricity e (such that 0 < e < 1) by 

a/q 2 — fr 2 
a 
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Kepler's second law 



* "for equal time intervals, a satellite will sweep out equal areas in its 
orbital plane, focused at the barycenter" 



Figure 2.2 Kepler's second law. The areas A } and A 2 swept 
ont in unit time arc equal. 

• Satellite takes longer to travel a given distance when it is farther away 
from earth (Explain ?) This can be useful ! 
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Kepler's third law 





* "the square of the periodic time of orbit is proportional to the cube 
of the mean distance between the two bodies" 



• Mathematically: 


o /i = 3.986005 x to 14 m 3 /s 2 is the earth's geocentric gravitational constant 
o mean distance is equal to the semi-major axis a 
o n is the mean motion of the satellite in radians per second 

• Orbital period P is related to n by P = 2zl 

• applies only to the ideal situation of a satellite orbiting a perfectly 
spherical earth of uniform mass, with no perturbing forces acting, 
such as atmospheric drag 
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Kepler's third law 



• Ex: Calculate the radius of a circular orbit for which the period is 1 
day. 



71 24x60x60 7.272 x 10 5 rad/s 


CL 


f 3.986005 xlO 14 1 § 
L (7.272X IO - 5 ) 2 J 


42, 241 km 
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Definitions of Terms for Earth-Orbitin 
Satellites 


Parameter 

Definition 

Sub-satellite path 

The path traced out on the earth's surface directly below the satellite 

Apogee 

The point farthest from earth (h a ) 

Perigee 

The point of closest approach to earth (h p ) 

Line of apsides 

The line joining the perigee and apogee through the center of the earth 

Ascending node 

The point where the orbit crosses the equatorial plane going from south to north 

Descending node 

The point where the orbit crosses the equatorial plane going from north to south 

Line of nodes 

The line joining the ascending and descending nodes through the center of the earth 

Inclination 

• The angle between the orbital plane and the earth's equatorial plane 

• measured at ascending node from the equator to the orbit, going from east to north 

• the greatest latitude, north or south, reached by sub-satellite path equal to 
inclination 

Prograde orbit 

An orbit in which the satellite moves in the same direction as the earth's rotation (0 -> 
90)/ (direct orbit)/ save launch energy 
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Definitions of Terms for Earth-Orbitin 
Satellites 





Parameter 

Definition 

Retrograde orbit 

An orbit in which the satellite moves in a direction counter to the earth's rotation (90 -> 

180) 

Argument of perigee 

CO 

• The angle from ascending node to perigee, measured in the orbital plane at the earth's 
center, in the direction of satellite motion 

Right ascension of the 
ascending node ft 

• Longitude of the ascending node is not fixed because the earth spins (Not a reference !) 

• Longitude and time are both required, but fixed reference in space is required also ! 

• Equinox (Vernal) -> when the sun crosses the equator going from south to north, and 
an imaginary line (line of aries) drawn from this equatorial crossing through the center 
of the sun points to the Equinox 

• The right ascension of the ascending node is then the angle measured eastward, in the 
equatorial plane, from this line to the ascending node 

Mean anomaly 

an average value of the angular position of the satellite with reference to the perigee 

True anomaly 

the angle from perigee to the satellite position, measured at the earth's center 
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Definitions of Terms for Earth-Orbitin 
Satellites * 
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Definitions of Terms for Earth-Orbiti 
Satellites 


Prcxjrade 

orbit 


Equator 



Figure 2.4 Prograde and retrograde orbits. 
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I = inclination 


A = Apoapais 


OB = 

Orbiting body 


Line of nodes 


CB = 

Central body 


Q = Argument 
of periapsis / 


Definitions of Terms for Earth-Orbiti 
Satellites 


Descending node 


Plane 

of reference 


P = Periapsis 


Ascending node 


»= Longitude 
of ascending node 


Vernal equinox 
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Definitions 

Satellites 


of Terms for Earth-Orbitins* 






Line of 
nodes 


Figure 2.5 The argument of perigee (0 and the right ascension 
of the ascending node Q. 
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Orbital Elements 



• Earth-orbiting artificial satellites are defined by six orbital elements 
referred to as the keplerian element set: 

o semimajor axis a and eccentricity e determines the shape of the orbit 

omean anomaly M 0 gives the position of the satellite in its orbit at a 
reference time (epoch) 

o Argument of perigee gives the rotation of the orbit's perigee point relative 
to the orbit's line of nodes in the earth's equatorial plane 

othe inclination and the right ascension of the ascending node relate the 
orbital plane's position to the earth 
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Two-line elements for NOAA-15 



dn 

year day J 

/ / \ 
i —u i j i 
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Orbital Elements 


• http://celestrak.com/NORAD/ele 
merits/ 

• Ex: from the following table, 
Calculate the semi-major axis for 
the satellite 


• Sol: 

Mean motion is given in rev/day 
as 


NN = 14.23304826 day 


-1 
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no = 2 x 7r x NN = 0.00104 s 1 

1 1/3 

= 7192.335 km 
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TABLE 2,1 Details from the NASA Bulletins (see Fig. 2.6 and App. 0) 


Line no. 

Columns 

Description 

1 

3-7 

Satellite number: 253 3 8 

1 

19-20 

Epoch year (last two digits of the year): 00 

1 

21-32 

Epoch day (day and fractional day of 
the year): 223.79688452 (this is 
discussed further in Sec. 2.9.2) 

1 

34-43 

First time denuatiue of the mean motion 
(rev/day 3 ): 0.00000307 

2 

9-10 

Jnc/iiifltion (degrees): 98.6328 

2 

13-25 

Rigfrt ascension of the ascending node 
(degrees): 251.5324 

2 

27-33 

Fccentrccify (leading decimal point assumed): 
0011501 

2 

3542 

Argument of perigee (degrees): 113.5534 

2 

44-51 

Mean anomaly (degrees): 246.6853 

2 

53-43 

Mean motion (rev/day): 14.23304826 

2 

6448 

ffeuotution num&er at epoch (rev): 11,663 
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Apogee and Perigee Heights 




• From the ellipse geometry, apogee and perigee radii are given as 
(derive !) 

• Apogee radius: r a = a(l + e) and Perigee radius: r p = a( 1 - e ) 

• Apogee and Perigee height are obtained by subtraction from earth radius 

• Ex: Calculate the apogee and perigee heights for the previous satellite 
parameters assuming mean earth radius is R = 6371 km 

• Sol: from the table e = 0.0011501 and from the previous ex. a = 7192.33 km 

r a = a(l + e) = 7192.33(1 + 0.0011501) = 7200.607 km 
r p = a{ 1 - e) = 7192.33(1 - 0.0011501) = 7184.063 km 

and the corresponding heights are 

h a = r a -R = 7200.607 - 6371 = 829.6 km 
h p = r p — R = 7184.063 — 6371 = 813.1 km 
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Orbit Perturbations 



Previous discussion assumes: uniform spherical earth, equilibrium of 
centrifugal force (satellite) and gravitational pull (earth) -> IDEAL ! 

Non-spherical earth (oblate spheroid) and atmospheric drag should be 
considered 


1. Non-spherical earth: Equatorial bulge and Poles flattening 

• The real mean motion 71 is modified to (in terms of n 0 , a, e and i): 

n = nnfl + K i( 1 ~ 1 - 5s ^ n r 2 ^ l A'i =66,063.1704 Km 2 

• u L ' a 2 ( 1 —e 2 ) 1 * 5 J 

• The orbital period (anomalistic period) is Pa = s 

• The mean motion is usually specified in the satellite parameters or obtained as the 
reciprocal of the anomalistic period in rad/s 

• Earth oblateness has negligible effect on semimajor axis a 
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Orbit Perturbations 



• Ex: A satellite is orbiting in the equatorial plane with a period from 
perigee to perigee of 12 h. Given that the eccentricity is 0.002, 
calculate the semi-major axis. The earth's equatorial radius is 
6378.1414 km 
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Solution Given data: e = 0.002; 2" = 0°; P = 12 h; K L = 660S3.1704 km"; a £ = 
6373.1414 km; fi = 3.936005 X 1014 m’W 

The mean motion is: 



2 IT 

n = — 

P 

= 1.454 X 10 s _1 



Assuming this-is- the same as n fh Kepler's third law gives 


■ - fcr 


= 26610 km 


Solving the root equation 


n — 


.£ 

s 


a L 


1 + 


K v ( 1 - 1.5sin 2 z) 


a 2 (l - tT) 


15 


= 0 


yields a value of 26,612 km 
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Orbit Perturbations 






• Regression of the nodes (due to oblateness of the earth) 

• the nodes appear to slide along the equator resulting is rotation of line of 
nodes and shift in the right ascension of the ascending node 

• If orbit is prograde -> nodes slide westward and if the orbit is retrograde 
nodes slide eastward (hence regression !) 

• rotation of apsides in the orbital plane 

• The mean motion n, the semimajor axis a and eccentricity e 
determines the amount of both effects through a factor K such that 

K = a 2 (i-e 2 ) 2 (rad/day) 

• The rate of change of ft with respect to time is ^ = -Kcosi 

fo/J/bg rate of change of u w i tim e is f = k (2 - 2.5 sin 2 i) 








Orbit Perturbations 






• At any time t, with respect to the epoch time t 0 , the angles are: 

(jj — ^{t — to) 17 = 17q + — to) 

• Thus, satellite drifts from elliptical path due to both effect (orbit itself 
is moving due to change in uj and n) 
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Orbit Perturbations 




• Ex: Determine the rate of regression of the nodes and the rate of 
rotation of the line of apsides for the satellite parameters specified in 
Table 2.1. The value for a obtained in previous Example may be used. 
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Orbit Perturbations 

Solution From Table 2.1 and Example 2.2 i - e = 0,0011501; MN = 

14.23304&26day -3 ; a = 7192.335 km. and die known constant if] = 66063.1704 km 2 



Converting n to rad/s: 


rt = 2ttNN 


From Lq. (2. II}: 


From Iiiq. C2.12}: 


K = 


ntC L 

<Al - eV 


= S_544 deg/day 


dil 

dt 


—JsT cos i 


= 0_9B1 deg/d^y 


From £iq_ (S_13): 


■Jlu 

dt 


= S — S _5 sin 2 e) 


= —2.904 d-eg/day 
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Orbit Perturbations 



• Ex: Calculate, for the satellite in previous example, the new values for 
oj and Q one period after epoch (use the table) 
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Orbit Perturbations 

Solution From Table 2. 1 : 

NN = 14.2,3304326 day- 1 ; fl* = 113.5534=; = 251.5324° 

The anomalistic period is* 



= 0.070259 day 

This is also the time difference since the satellite has completed one rev¬ 

olution from perigee to perigee^ Hence: 

a = tit, + - t 

dt 

= 251.5324 -I- 0.931(0.070259) 

= 251.601° 


iii = ii > 0 +- 




dt 

= 113.5534 
= 113.349° 


t - o 

+ ( — 2.903) (0.070259) 



10/10/2016 


Copyright © 2016 - Mahmoud Selim 


25 








Orbit Perturbations 






• Equatorial ellipticity (of the order of 10 -5 ) (because earth is not 
perfectly circular in equatorial plane) results in gravity gradient that 
affects satellites in geostationary orbits to drift rather than being fixed 
to earth 

• This effect is handled by station-keeping maneuvers 
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Orbit Perturbations 



2. Atmospheric drag 

• Affects near earth satellites (below about 1000 km) 

• The effect is highest at perigee reduction in velocity apogee 
height decrease successively -> e and a decrease 

• Significant for long time intervals 


Approximate expression for majoraxis change at time t 


a = a o[ 


n 0 


no+n 0 (t-to) 


] 2/3 


The change in mean anomaly is 

SM = #(t-t 0 ) 2 
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Calendars 



• Tropical year contains 365.2422 mean 
solar days (relative to mean sun) is 
divided into 365 days 

• The extra 0.2422 day is corrected by 
the "leap year" Julian Ceaser -> 
Julian Calendar 

• Gregorian Calendar is in use today -> 
Pope Gregory XIII by solving a problem 
by missing out 3 days every 400 years 

• Sidereal day: one complete rotation of 
the earth relative to the fixed stars (1 
mean solar day = 1.0027379093 mean 
sidereal day) 
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Universal time 



• UTC Universal Time Coordinated 

• Reference time for setting clocks 

• The fundamental unit is the mean solar day (24h,60m,60s) 

• Epoch time is given in terms of UTC 

• Similar to Greenwich mean time (GMT) and Zulu (Z) time 

• UT is expressed in two forms: as a fraction of a day and in degrees 


UT d 


ay 


_ 1 f hm,r*Q I minutes . seconds^ 

94 ['vULLf O “1 "I 


60 


3600 J 


UT° = 360° x UT day 
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The orbital plane 


• From geometry, magnitude of 

vector r is 

r = q( i-e 2 ) 

1 -\-ecosu 

• How to find v (true anomaly) 
which is varying with time ?! 

• First, find the mean anomaly M 
at time t as 

M — n(t — T p ) 

• But, Tp = to - M* 

• Thus, M = Mq + n(t — to) 
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The orbital plane 



Secondly, solve Kepler's 
equation for an intermediate 
variable E (Eccentric anomaly) as 

M — (E — esinE ) = 0 

Solved using iterative method ! 
Find v using Gauss's equation: 

tan | = yj jz^tan^ 

The magnitude of vector r can 
also be obtained from E as 

r = a(l — ecosE) 



Figure 2.8 Perifocal coordinate system (PQW frame). 
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The orbital plane 



• For near-circular orbits, v can be 
obtained directly from M as 

v = M + 2 esinM + |e 2 sin2M 


Vu. 
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The orbital plane 




• Ex: For a given satellite, e = 9.5981 x 10 3 _and mean anomaly is 
204.9779 deg v the mean motion is 14.2171404 rev/day, a = 7194.9 
Km 

Calculate the true anomaly, the magnitude of the radius vector 5 s 
after epoch 


10/10/2016 


Copyright © 2016 - Mahmoud Selim 


33 



The orbital plane 

Solution The rotation in radians per second is 



14.2171404 X 2w 
n =- 

S6400 

= 0.001034 rad/s 



The mean anomaly of 204.9779 a h in radians is 3.57754, and 5 s after epoch the 
mean anomaly becomes 

M = 3.57754 + 0.001034 X 5 

= 3.5827 rad 

p= 3.5827 + 2 X 9.5981 X 10 -3 X sin 3.5827 

+ - X (9.5981 X 10 -3 ) 2 X sin (2 X 3.5827) 

4 

= 3.5746 rad <= 204.81°) 

Applying Eq. (2.23) gives r as 

7194.9 X<1 -9.5981 s ) X 10 6 

r -'- 

1 + 9.5981 X 10" 3 X cos 204.81 

= 7257.5 km 
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